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CHARTS GIVING CRITICAL COMPRESSIVE STRESS OF CONTTNUOUS 
FLAT SHEET DIVIDED INTO PARALLEL.OGRAM- SHAPED PANELS 
By Roger A. Anderson 


SUMMARY 


Charts giving the congressive-hnckling-stress coefficients for 
sheet panels of a shape occurring in swept-wing plan forms are presented. 
The panels are assumed to he a part of a continuous flat sheet divided 
hy nondeflecting supports into parallelogram- shaped areas. The stability 
analysis was performed hy the energy method and the results show that, 
over a wide range of panel aspect ratio, such panels are decidedly more 
stable than equivalent rectangular panels of the same area. 


INTRODUCTION 


A hl^ly desirable characteristic for hl^- speed fll^t is to have 
aircraft outer surfaces that remain free of buckles or waves under all 
normal fli^t conditions. Whether these surfaces remain smooth under 
fll^t loads is determined by the skin thickness and the arrangement 
and rigidity of the internal supporting structure to which the skin is 
attached. In the past, the supporting structure generally divided the 
outer skin into an array of approximately rectangular panels, but 
present-day swept- and delta-wing plan forms call attention to the fact 
that skin panels of other shapes such as parallelograms- may occur. 

The present paper considers the stability under ccm 5 )res&ive stress 
of continuous flat sheet divided by nondeflecting supports into an array 
of paiallelogram- shaped panels. Wide ranges of panel skewness and aspect 
ratio were investigated, eind two orientations of the parallelogram- shaped 
panels with respect to the direction of the applied stress were con- 
sidered. The results of the analysis are presented in the form of charts 
of theoretical buckling-stress coefficients as a function of panel skew- 
ness and aspect ratio. 
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SYMBOLS 


mutually pei^iendlcular directions, rectangular 
coordinate system 

direction parallel to x-direction, skew coordinate system 

arbitrary direction, skew coordinate system 

skew angle between y- and y* -directions, measured 
positive clockwise from y-directlon, degrees 

panel dimension In x- or x’ -direction 

panel dimension in y-direction 

panel dimension in y’-dlrection 

panel aspect ratio (a/b) 

plate thickness 

deflection normal to plane of plate 
plate flexural stiffness 
Poisson's ratio for plate material 
general notation for stress 
stress acting in x-directlon, 
stress acting in y-dlrectlon, 

critical-stress coefficients 


con5>ression positive 
con5>reBsion positive 


( (rtb^. ^ 



internal bending energy 


T 


external work of app 3 J.ed stress 
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j,m,n Integers 

amn ,~bnm Fourier coefficients 


DESCRIPTION OF PROBLEM 


An idealized arrangement of the supporting structure possible for 
svept and delta wings is shown in figure 1^ wherein a part of a con- 
tinuous flat sheet infinite in length and width is shown divided by non- 
deflecting supports into an array of parallelogram- shaped panels, of 
which rectangular panels are a special case. For the case of rectangu- 
lar panels the assun5)tion is usually made that, if the supporting mem- 
bers are torsionally weak, the plate bending moments arising during 
plate buckling under edge con5)ression stress are negligible at the panel 
boundaries and each panel therefore may be treated as an isolated rec- 
tangular plate with sin^jly supported edges. Stability data for such 
panels are available in standard textbooks, such as reference 1. 

When a simply supported rectangular plate in con^yresslon is skewed 
into the shape of a parallelogram, a certain increase in stability due 
to the shape change is achieved,' as shown by the numerical exanqples in 
reference 2. A further increase in stability is possible when the 
skewed panels are part of a continuous sheet because of the restraint 
the adjacent sheet panels lng)ose on each other during the formation of 
a continuous buckle pattern. The effect of this type of restraint is 
illustrated in reference 3 foi* 'tlie case of continuous sheet divided into 
square panels subject to shear stress. Restraint due to sheet con- 
tinuity would be present in the case of parallelogram- shaped panels 
regardless of whether the loading is edge cong)resslon or shear. 

Two orientations of the parallelogram- shaped panels with xespect 
to the principal direction of the con5)ressive stress due to wing bending 
are considered in the present paper. The two loading conditions are: 
stress acting parallel to a set of panel sides and stress acting perpen- 
dicular to a set of panel sides. Both loading conditions are shown in 
figure 1. ' 

In the analysis the assung)tlon is nade that the supporting members 
are rigid enough to prevent deflection of the sheet at the panel 
boundaries but offer no restraint to rotation. A quantitative analysis 
of the stiffnesses required of supporting members (ribs, stiffeners, 
shear webs, etc.) to prevent deflection around the edges of parallelogram- 
shaped sheet panels is beyond the scope of this paper, but indications 
are that they would be somewhat higher them the stiffness required to 
support equivalent rectangular sheet panels. 
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In order to establish a staMllty criterion for the configuration 
analyzed, the energy method of analysis is used, and an approximation is 
made for the deflection of the plate expressed in skew coordinates (see 
appendix A) . The deflection function used leads to the exact solution 
of the differential equation of equilibrium of the plate loaded in com- 
pression for the special case of rectangular panels and also for the 
special case of equal-sided panels of arbitrary skewness subjected to 
conqjressive loading perpendicular to one set of sides. The accuracy 
of the rest of the data is indicated, "sdiere feasible, by conq)arlson 
with the results of more accurate energy solutions, idiich are given in 
appendix B. 

RESULTS AND DISCUSSION 


The critical conq)ressive stress for parallelogram- shaped plates may 
be given by the formula used for rectangular plates 

(Tt = k 

b 


where the dimension b is the peipendlcular distance between the sup- 
ports alined in the x-dlrection, as in figure 1. The critical-stress 
coefficient k depends on the panel aspect ratio p (defined as a/b) , 
the direction of the applied compressive stress, and the magnitude of 
the skew angle 0. 

For loading in the x-dlrection (stress acting parallel to a set of 
sides), the chart for k^ is given in figure 2. Rote should be taken 
that, for a given value of p, the buckling coefficients indicated for 
each curve are associated with panels of equal area. In order to facili- 
tate association of the curves of figure 2 with the geometry of the 
panels, figure 3 bas been prepared. In this figure, sketches of the 
panels of aspect ratios of 1, 2, and 3 a-i’s presented along with four of 
the curves of figure 2. In each sketch the mode of buckling for the 
panel is indicated. Figures 2 and 3 indicate that the stability of the 
skewed panels is definitely increased relative to equivalent rectangular 
panels of the same area over a wide range of panel aspect ratio. This 
increase in stability is due partly to the change in panel shape but is 
caused mainly by the restraint lii 5 »osed on the mode of buckling due to 
the presence of adjacent skewed peinels. At aspect ratios of 4.5 and 
greater, this restraint has largely disappeared and the coefficients 
approach the value 4. 

The curves presented in figures 2 and 3 were derived by the energy 
method of analysis with the use of an assumed deflection function capable 
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of giving an idealized representation of the panel ‘buckling modes. The 
buckling coefficients thus obtained represent upper limits to the true 
coefficients. In order to evaluate the unconservativeness of the curves 
presented in figures 2 and 3^ check calculations were made by energy- 
solutions of greater accuracy. The results of these calculations, -which 
are believed to be less than 1 percent unconservative, are indicated by 
the circles near the curves in figures 2 and 3* Because of the com- 
plexity of the buckle patterns, check calculations over the entire range 
of aspect ratio were not feasible to make, but it seems reasonable that 
the remaining parts of the curves would show about -the same degree of 
unconservativeness as in the regions where check calculations -were made. 
The numerical values of the circles in figures 2 and 3 are listed in 
table 1. 

For loading in the y-direction (stress acting perpendicular to a 
set of sides), the chart for ky is presented in figure These 
curves show that a single buckle pattern (characterized by an alternate 
in and out buckling from panel to panel) exists throu^out the range of 
aspect ratio investigated except for the expected change in mode at 
aspect ratio 1//2 for the 0=0° case. Again, the effect of con- 
tinuity bet-ween psuiels on -fche buckling coefficient has largely disap- 
peared at aspect ratios of and greater. The sq.uares on these curves 
indicate the points at -vdiich the derived curves pass throu^ -the exact 
value for the buckling coefficient, and -the results of energy solutions 
of hl^ accuracy indicated by the circles are in excellent agreement 
with the remaining parts of the derived curves. (See -table 1 for numeri- 
cal values. ) 

Under the simultaneous action of con^iressive stresses in the x- and 
y-directions, the critical- stress combinations for a given peinel can be 
shown by means of an interaction curve. A family of such curves is pre- 
sented in figure 5 for an array of equal- sided panels for -various angles 
of ske-wness of -the panels. It is rather Interesting that for equal- 
sided panels the s-tability criterion for a Uy type of stress acting 
alone is Independent of 0 and gives a constant -value of 4 (exact) for 
the buckling coefficient. 

The interaction curves for 0 = 45° and. 60° in figure 5(a) have 
been adjusted in -the regions -vdiere -the more accurate energy solutions 
showed -that the energy solution using an idealized represen-tatlon of 
the buckling modes -was unconservative. The relation be-tween -the more 
accurate check calculations and the approximate solution is shown in 
figure 5(t)* The discontinuities in slope of -the curves for 0 = 45° 
and 60° in -this figure indicate that several changes in buckling mode 
take place in passing from a cfy type of stress to a "type of stress. 

The smooth curves in figure 5(a), however, were ob-tained by singly 
fairing a curve -throu^ -the check points of figure 5(b) . The numerical 
values of the check points are included in table 1. 
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CONCLUDING BEMASKS 


The charts presented show that the stahillty of flat continuous 
sheet divided hy supporting memhers into skewed panels is definitely 
increased in relation to the stability of equivalent rectangular panels 
of equal area over a wide range of panel aspect ratio. This increase 
in stability may be attributed partly to the nonrectangular shape of 
the panels but is caused mainly by the restraint ingposed on the mode of 
buckling due to the presence of adjacent skewed panels. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., April 5> 1951 
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APPENDIX A 

APPROXIMATE SOLUTION FOR CRITICAL COMPRESSIVE STRESS OF 
FLAT PLATE CONTINUOUS OVER SKEWED SUPPORTS 


Consideration of a flat plate supported and loaded as shown in 
figure 1 reveals that, for many combinations of plate aspect ratio and 
skew' angle, a rather complex deflection pattern. Tuight be expected when 
buckling occurs. The associated plate differential equation of equi- 
librium written in skew cooid.inates (x’,y*) is 


Bx’^ 


+ 2(1 + 


2 sin20) 


S^w 



k sin 



5 ^w \ 
3x*dy* V 


ii cosV - 
° d. 



Sx* Sy’ 


sin 



This equation further Indicates that for arbitrary panel dimensions an 
exact solution for the deflection w is not likely to be found. It is 
interesting, however, that for the particular case of equal-sided skew 
panels under congjression stress acting peipendicular to one set of 
sides, the familiar deflection function ' 

w = sin sin 

a a 


written in skew coordinates is the exact solution to the differential 
equation. The resulting stability criterldn is Independent of the skew 
angle ^ and the buckling- stress coefficient ky is equal to k, or 


0 yt = 


k«^D 

b^ 


The points where this solution applies have been noted on the curvea of 
figure k. In order to calculate the buckling-stress coefficients for 
other panel dimensions, the energy method can be used -with an assumed 
expression for the plate deflection. 
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Choice of deflection function .- The foregoing exaii5)le Indicates 
that, for skewed panels con^iressed In a direction perpendicular to one 
set of sides (y-dlrectlon In fig. 1 ), a single skewed sinusoidal mve 
In each panel Is a good assumption for a deflection function for panels 
of eq,ual or nearly equal sides. This assumption Is also made for panels 
of larger aspect ratio. For conqiresslon loading parallel to one set of 
sides (x-dlrectlon In fig. 1 ), however. It Is evident that for most of 
the range of aspect ratio, several buckles will form In each panel. 

Each panel therefore won Id contain contours of zero deflection, or nodal 
lines, between buckles. Since continuity of deflection and slope must 
be preserved at the peinel boundaries, the nodal lines must run continu- 
ously across panel boundaries. For example, skewed panels of large 
aspect ratio under stress In the x-dlrectlon would buckle In such a 
manner that a system of rather uniformly spaced nodal lines would 
develop between buckles and would run continuously In the general direc- 
tion of the y- coordinate. This characteristic behavior. Idealized In 
the following sketch. 


y y' 



a 

would be expected for con^jresslon In the x-dlrectlon whatever the skew 
angle or aspect ratio. Dashed lines, such as a-a, represent nodal lines. 
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In the sketchy the buckle pattern is represented as being sym- 
metrical about a panel midpoint and repeating in every panel. However, 
in general, a repetitive buckle pattern can be assumed to take place 
in a group of panels of unknown dimensions Ha in the x-direction and 
Hb’ in the y’-directlon, idiere H and M are integers. 

A representation of such a pattern can be built up in two steps. 
First, zero deflection along all panel boundaries can be provided by 
the function ^ 


^1 


sin 


pitx’ 

Ha 


sin 


Mb* 


(Al) 


Tdiere the ratios p/H and q/M are integers. For stress in the 

y-direction, equation (Al) with §■ = 1 ^ = 1, 2, 3> • • • is a 

suitable approximation for the buckle pattern which occurs over a wide 

range of panel aspect ratio and skew angle 0. Similarly, the buckle 

pattern occurring in sli^tly skewed panels of moderate aspect ratio 

under stress in the x-directlon is represented by eq.uatlon (Al) with 

^ = 1, 2, 3^ • • • and — = 1. For this special case, the nodal 
H M 

lines a-a in the sketch are parallel to the y’ -direction. In order to 
represent all other orientations of the nodal lines a-a and to permit 
an arbitrary uniform spacing, W]_ can be multiplied by some function 
such as V2f vhere 


W£ = cos 




(A2) 


The lines along which W 2 = 0 have a slope determined by the 
ratio Ha/Mb* and a spacing determined by the ratio J/H, where j 
is an Integer and may taJce on all values 0, 1, 2, 3^ • • . . 

By multiplying wq_ by W 2 , the desired deflection function 


w = sin 


pnx* 

Ha 


sin 


<lrty* 

Mb* 




(A3) 


is obtained. This function is capable of satisfying the idealized 
physical characteristics of the buckle pattern for einy panel aspect 
ratio and skew angle upon proper choice of the integers p, q, H, M, 
6ind j. Actually, it is necessary to consider only the plus sign in 
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equation (A3) under con^reseion loading since the plus- sign case descrihes 
all natural orientations of nodal lines. Another form of equation (A3)^ 
more convenient for calculation purposes, is 


w = 


sin 


p«x* 

Na 


. qity* / Jrtx' 

sin - ■ ■ ■(cos 

Mb‘ \ Na 


Jjty’ 

cos 

Mb* 


sin 


Jjtx’ 

Ka 


sin 


Mb*/ 


(A4) 


Energy expression .- The internal energy V of the plate anfl the 
external vork T done hy the loads, eapressed in rectangular coordinates, 
are 



Conversion of the energy expression into the skew coordinate system 
(x*,y*)' is acconpllshed with the transformation equations 


X* = X - y tan 


y* = 


y 

cos ^ 


where the positive directions of the coordinates are shown in the 
following sketch: 
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As shown in equations (A 7 ) and (A 8 ), the Integrations are performed 
over a repetitive huckle pattern. 


Stahlllty criteria .- By substituting equation (A 4 ) for w into 
equations (A 7 ) and (ASJ and hy setting V equal to T, the following 
stability criterion is obtained: 



and 


A 


2 
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A somewiiat more convenient form of equation (A9) for computational 
puii)OBes Is ol)talned If the skewsd-plate aspect ratio Is defined as 
a/b to conform with the usual notation for a rectangular plate of equal 


area, 
kx’ cos 


"h B ’ 

If the substitution b’ is made, then — - — 3 - 

cos p cos p 

^ = kx, ky’ cos^ = ky, and equation (A9) becomes 






The critical values of the buckling- stress coefficients k^j and ky 
are determined by the values of the ratios p/N, 9/^^ J/M that 

make the coefficients a minimum. 

If which is equivalent to using equation (Al) for a 

deflection function, equation (AlO) reduces to 

, e I J , 4 g I , D (AU) 


This equation is the exact stability criterion for a rectangular plate 
in compression when ^ is set equal to zero. The curves of buckling 
coefficient plotted against aspect ratio for stress in the y-dlrection 
in figure 4 were conq)uted from equation (All) by setting kx = 0. Over 
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th.e range of aspect ratio covered, the lowest huckllng coefficients were 
associated with the parameters ^ ^ except for the part of the 

curve for 0 = 0° in the region i < P < for which the huckle pat- 


tern is defined hy — = 1, — = 2. 

N M 






For the 

and for £ = 
W 


special case 

— = 1 arifl k-v 


P = -T, Tdiich specifies equal-sided pemels, 

cos p 

= 0, equation (All) reduces to 


ij-(l + sin^) = ky(l + sin^) 


or 



This result is the one obtained previously hy substituting 
V = sin — — sin — ^ into the plate differential equation. 

A determination of the huckling-stress coefficients for stress in 
the x-direction from equation (AlO) may he made when a huckle pattern 
has heen specified hy the iJarameters p/N, j/N, q/M, and J/M. For 
a given panel shape, several huckle patterns must he Investigated to 
determine the pattern associated with the lowest huckllng-stress coef- 
ficient. The adnimum con^juted coefficients thus obtained are plotted 
in figure 2, and the parameters p/N, j/N, q/M, and J/M defining 
the huckle pattern in each part of these curves are tabulated in table 2 . 
From the table, the data for p = and P = 3*2 were used to sketch 

the huckle pattern shown in a previous section "Choice of deflection 
function." Similarly, the huckle pattern associated with any other 
panel configuration covered in figure 2 can he sketched hy substituting 
the appropriate parameters from table 2 into the general deflection 
function, equation (A 3 ) . This procedure was followed in preparing 
figure 3 . 

The Interaction curves in figure 5 were derived hy trial-and-error 
solution of equation (AlO) hy the process explained in the previous 
paragraph. The parajneters defining the huckle pattern in each segment 
of these curves are tabulated in table 3 * 
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APPENDIX B 

MOEE ACCUMTE SOLUTION FOR CRITICAL COMPRESSIVE STRESS OF 
FLAT PLATE CONTINUOUS OVER SKEWED SUPPORTS 


The results ohtained in appendix A serve as a useful guide in 
making a more accurate energy solution for the critical stresses 
because the basic characteristics of the mode of buckling for a number 
of skewed panel configurations have been determined. These data are 
given in tables 2 and 3- A study reveals that a number of the buckle 
patterns can be classified rather simply as being either symmetrical 
or antisymmetrical about the midpoint of each panel and repeating in 
each panel or over two panels. The solutions in this appendix are 
confined to these cases. Representation of such patterns is feasible 
with, two-dimensional sine and cosine series and permits an energy 
solution for the critical stresses to an arbitrarily hi^ degree of 
accuracy. The critical stresses thus obtained provide a check on the 
results in appendix A. 

Deflection functions .- The trigonometric series representing the 
deflection were chosen to satisfy term by term the required conMtlons 
of zero deflection around the panel boundaries and continuity of slope 
between panels. The four buckling configurations investigated are as 
follows : 

Symmetric buckling, periodic over a, b’, 

w= T X sin 2fl sin + 

m=2, k, 6 n=2, 4, 6 ^ 


m=l,3,5 n=l,3,5 


cc 


bjjm I cos 


(m - l)rtx* 


cos 


(m + l)mx^ (n - l)jty» _ (n + l)jty» 

a ii b» b» _ 


(Bla) 
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Symmetric ‘buckling, periodic over 2a, 2b*, 


00 00 


V = 


. mjtx’ . nity’ 
ajuQ sin — r — sin — — + 


m=l,3,5 n=l,3,5 


00 00 


m=2,k,6 n=2,4, 6 


, , (m - l)jtx* 

^mn|cos 


cos 


(m + l)itx 


E 

— cos — 


- l)ity* (n + l)rty* 

— cos • 


"b* 


Antisymmetric 'buc kli ng, periodic over a, 2b*, 


00 00 


. mjtx* . nirv* 

V = ^ ^ amn sin sin — ^ 1 - 

m=2,k,6 n=l,3,5 ®' "b* 


00 00 


r“ 


(m - l)rtx* 

'O — - — 

m=l,3,5 n=2,l^,6 


■bjnji cos 


cos 


[cob - cos :LA)jgl| 

a JL J 


Antisymmetric 'buckling, periodic over 2a, ‘b*, 

00 00 

'ST X . nuDc’ . njty* , 

V = ^ ^ a^BixL Bln + 

m=l,3,5 n=2,4,6 ^ 

m=2,it,6 n=l,3,5 L 

[cob (a - i)*v» _ cos Ln + _i)jar»7| 

JL J 


(Bl'b) 


(Blc) 


cos 


(Bid) 



u 


Substitution, of these deflection functloais Into the energy expresslonB (equations (A7) and 
(a 8)) and minimization of the potential-energy function (Y - T) vlth respect to amn and hnin 
leads to the desired stability eq,uations. 

Symmetric buckling, periodic over a, b* .- For symmetric buckling, periodic over a, b‘, 

■ • « •!. .1 n ^ ^ ^ ^<1 trj-, ^>1 /'Q^Q^ -f O 

■cue lHI3JlX1>e BCC U1 B'OtiU±M±OJ CliUCtUJ-UJ-LB U.O±XV^=^J. AUJ^V^ uj-wj-k y^j^-L^y U-C 

n/" m2 2 P ^ fl< 

0 - amn^ + Wi jb(nH-l) (n+l) ' D(m+l)(n-l) " D(m:i)(n+l) + ^(m^l) (n-ljj o . | k^l tP^n «! 


(m »■ 2, i*-, 6, . . .) (n 2f k, 6, . . .) (B2a) 


. _ ,L/ .s\/. f m + 1 \ ^ ^ \ 2 & .>_^j ^ L /_ . t\/— ^\_ 




(HU-1) (d-1) 




•,\2 e 




l^(ln - 11 (a+ l)a^~ [ I“ ~ r*-) - +6(n + l)^-E t. 

' cta*^ c ^ 


tun + V(m -l) (n- l)a/„ i\/„ ni • + B^n- 1)^-2- k- 

' coi’.^ 


tan i—Zb^f^ r ,\ 1 ^ + 

'Lpi^coo*^ 


p2(ii - 1)*^+ 2(1 + 3 ton^) (n - l)^(-2 + 1) - (k 3 + ky tan2|i)(«2 + l) . p^(n - 1)^ - a^nn-a) + D^ + 2(1 + 3 t«i^) (n + 1)2 (b 2+ i) _ 

(kx+ky tan2jj)(]a2 + i) -p%£y(n + l)^ - 2b(^g)„ 3^(0*^ + 6ti^ + 1) + 2(1 + 3 taa^) (m - l)2(n2+ 1) - (k^+ky taa^)(m -1)^- 3^(n2 + 1^ _ 

a>(a.E)n|^^^+ + &2+ 1) + 2(1 + 3 tan^) (■ + 1) 2(n2 + 1) . (i^ + taa^) (b+ l) ^ - p=ky(n2+ 

2(1 + 3 tan^)(m-l)®(n-l)^- (ki + fcy tan^)(n -1)®- p^n- l)J + b(jn. 2 )(iixg) -i^t^ + p2(n +l)^ + 2 (l + 3 tan^^) (■+ l)®(n + 1)® - 

“n I— . 

(kx+ky t*m2ji)(*i+l)^-p^y(n + l)^ +t(,^2){j^^)|^^^+p2(n + l)* + £(l + 3 tan^)(»-l)^n + l)^- (kx+l^j tan^^) (a -1)2 -p®ky(n+l)^ + 

’’(»f2)(D-2) tP^(p - !■)**'+ g(H- 3 tan®^)(a + l)^(n-l}®-(kjc+ky tim2fj)(« + i)2 -p^(n- 1)^ 


(■ - 3, ?, T, • . .)(= - 3. 5, 7» • • •) (B2b) 
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0 = binBin - 8(n - l)a2(n-l) 


P 'cos^ 


p(n - 1)^ - I ky 


tan ^ + 


8(n + l)a2(n+i) 


ii— + p(n + 1)^ - I ky 

p cos^ 




16 

P^cos^0 


3P^(n +1) + 8(1 + 3 tan^)(n + l)^ - ^(kjj + ky tan^) - 


3 P%(n + 1 )^ 


" '*^-l(n-2) 


16 ^„2, ^.4 

'+ 3 P (a - 1) + 


p^cosV 


8(1 + 3 tan20)(n - l)^ - ^(kx + ky tan^) - 3P^ky(n - l)' 


2b 


3n 


16 


P^cos^0 


+ P^(n^ + 6n^ + 1) + 8(1 + 3 tan^)(n^ + 1) - 


4 (kx + ky tan^) - p^ky(n^ + 1) 


+ ^3(n+2) 


16 „2/ ,\4 

+ p (n + 1)^ + 


P^cos^0 


8(1 + 3 tan^)(n + l)^ - 4 (kx + ky tan^) - p^y(n + 1)' 


^3(n-2) 


16 


— j— + p2(n - 1)^ + 8(1 + 3 tan^0)(n - l)' 

P^cos^0 


4(kx + ky t€in^) - p^ky(n - 1)‘ 


(n = 3 , 5 , 7 , • • • ) (B2c) 
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0 = VLBml - - l)a(m.i)2 


- 1) + 4b - i k 
|l.cos2^ 2 y 


tan 0 + 


8(m + 


+ 1,^ - I ky 

P coa^ 




3(m - ir ^ 
P^cos^0 


l6p^ + 8(1 + 3 tan^)(m - l)^ - 3(kx + ky tan^)(m - l)^ - 4p^ky 


^(m+2)l 


+ 8(1 + 3 tan20)(m + l)^ - 


P^cos^0 


3(kx + ky tan^)(m + l)^ - ^P^ky 


- 2bni3 


+ 6m^ + 1 , 2 


1 p^cos^0 


+ l6p^ + 


8(1 + 3 tan^)(m2 + l) - (k^ + ky tan20)(m^ + l) - 


Xn^2)3 




^ + l6p^ + 8(1 + 3 tan^) (m + 1)^ - 


P^cos^0 


(kx + ky tan20)(m + l)^ - ^P^ky 


^(m-2)3 


(a - + i6p2 + 


P^cos^ 


8(1 + 3 taji^)(m - 1)^ _ (kx + ky tan^)(m - l)^ - 4 p^- 


-] 


(m = 3, 5, 7, . . .) (B2d) 



20 


MCA TN 2392 


0 — ~b -| 1 B -| 1 + l6& pp / 


COB 




+ ^ tan 0 - 


4 2 

+ 123 + 


3 ^cos ^0 


8(1 + 3 tan^) - (kx + % tan^) - 33^^^^ 


- 4 l 3 


31 


12 , „2 

+ + 


3 2 c 08^0 


8(1 + 3 tan^) - 3 (t;x + ky tan^) - 3^1^ 


+ irt) 


33 


3 ^cos ^0 


+ 43 ^^ + 


8(1 + 3 tan^) - (kx + ky tan^) - 3^- 


(B2e) 


where 


Athti — 




3 ^cos^ 


+ n^3^ + 2m2n2(i + 3 tan20) - m2(kx + ky tan^) - n232ky 


Bum. 


. ^ + &2 + 1 ) + 

p2cosV 

8(1 + 3 tan^)(m2 + l) (n^ + l) - 


4(kx + ky tan^) (m^ + l) - 43^ky(n^ + l) 


Bln = — — + 63 ^(n^ + 6n2 + l).+ l6(l + 3 tan^) (n2 + 1) - 

3^008 V 

8(kx + ky ten^) - 63 ^y(n^ + 1) 
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®ml 


= 6(m + 6m + l) ^ ^ ^ 2 tan^)(m2 + l) - 


p2cos^0 


6(kx + ky tan^) (m^ + 1) - 8^^ky 


48 

Bii = j— + ltB^^ + 32(1 + 3 tan 20 ) - 12 (kx + % tan^^) _ 12p2t 


^^cos ^0 


The stability criterion is obtained by equating to zero the determinant 
formed by the coefficients of amn and brim in equations (B2) . The 
numerical results calculated from this set of equations may be found in 
table 1. In each case four of each of the coefficients amn and bnmj 

vhich gave an ei^th- order determinant, were used in the calculations to 
insure an adequate representation of the buckle pattern. 

Symmetric buckling, periodic over 2a, 2b* .- For symmetric 

buckling, periodic over 2 a and 2 b’, the infinite set of stability 
equations derived from the function (Bib) is the same as that derived 
from the function (Bla) except for a change in the values of m and n 
involved. For the buckle pattern now under consideration, equations (B2a) 
exist for m = 1, 3, 5^ • • • an<i n = 1, 3^ 5^ • • • • Similarly, equa- 
tions (B2b) exist for m = 2, If, . . . and n = 2,lf, . . ., and equa- 
tions (B2 c), (B2d), and, (B2e) do not exist. Note should be taken that 
coefficients with the subscripts m - 1 and n - 1 in equations (B2a) 
drop out for m = n = 1 since coefficients with a zero subscript do not 
appear in the deflection function (Bib) . For the same reason coefficients 
with the subscripts m - 2 and n - 2 should be dropped from equa- 
tions (B2b) for m = n = 2. 

The numerical results calculated for certain panel configurations, 
in -vdiich a symmetrical buckle pattern repeating over 2 a and 2 b’ is 
associated with the lowest buckling loeid, are given in table 1. For the 
panels in -vdiich 0 = 30 ° and 45°, four equations involving a-]-] , 

a^q, and bpp were employed in the calculations. For the panels of 
60 ° skew, four of each of the coefficients a mri and httm were required 
to provide an adequate representation of the buckle pattern. 

Antisymmetric buckling, periodic over a, 2b’ .- For antisymmetric 
buckling, periodic over a and 2 b*, the infinite set of stability equa- 
tions derived from the function (Blc) may be written directly from 
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equations (B2a) to (B2c). For this huckle pattern, equations (B2a) 
exist for m = 2, 4, 6 , . . . and n = 1, 3^ 5^ • • •; equations (B2b) 

exist for m = 3 ^ 5^ 7^ • • • ani n = 2, Ij-, 6 , . . . ; equations (B2c) 

exist for n = 2, k, 6, . . and equations (B2d) and (B2e) do not exist. 

The buckling coefficient k^ for loading in the x-direction was calcu- 
lated for the panel configuration 0 = 60 °, P = 3 for which this huckle 
pattern applies, and a value k^ = 6*95 was obtained. This result was 
obtained when four a’s and four b*s were used in the deflection 
series (Blc) . 

Antisymmetric buckling, periodic over 2a, b* .- For antisymmetric 

buckling, periodic over 2a and b*, the infinite set of stability 
equations derived from the function (Bid) may be written directly from 
equations (B2a) to (B2d) . For this buckle pattern, equations (B2a) 
exist for m = 1, 3 , 5^ • • • and n = 2, k, 6 , . . .; equations (B2b) 

exist for m = 2, k, 6 , . . . and n = 3^ 5^ If • • • ; equations (B2c) 

do not existj equations (B2d) exist for m = 2, 6 , . . . ; and equa- 

tion (B2e) does not exist. 

Calculations were carried out for p =. 60°, ^ = 1.2, a panel 

configuration for which this buckle pattern applies. An adequate 
representation of the deflection was obtained when four a*s and 
four b*s were used in the deflection series (Bid) which gave a value 

of kx = 16 . 22 . 

As was pointed out in the beginning of this appendix, these more 
accurate solutions were carried out only for those panels for which the 
mode of buckling was relatively sin^jle and clearly Indicated by the 
results of appendix A. In the absence of such information it would be 
necessary to investigate all conceivable modes of buckling for a given 
panel configuration. 
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TABLE 1.- CEITICAL COMBIHATiaNS OF CTx AND ffy TYPES OF STRESS ACTING 
ON CONTINaoaS FLAT DIVIDED INTO PAEALLELOGEAM- SHAPED’ PANELS 

(DATA FOE CIRCLES IN FIGS. 2 TO 5) 

^ b2 ^ 



p 

Stress combination 

Buckle pattern 


% 

Symmetry 

(1) 

Periodici-ty 



0 = 

30° 



0.5 

0 

19.33 

S 

2a 

2b* 

.577 

9.60 

0 

S 

’2a 

2b* 

1.0 

6.7k 

0 

s 

2a 

2b* 

' 1.155 

6.62 

0 

s 

2a 

2b* 

1.155 

0 

4 

s 

2a 

2b* 

3.0 

0 

i. 4 i 

s 

2a 

2b* 



0 = 

45 ° 



0.5 

24.12 

0 

s 

2a 

2b* 

. 6 

0 

14.32 

s 

2a 

2b* 

1.0 

0 

6.499 

s 

2a 

2b* 

1.0 

11.46 

0 

s 

a 

BOH 

i. 4 iij- 

0 

4 

s 

2a 


i. 4 i 4 

5.99 

2 

s 

2a 


i.iaii- 

7.46 

1 

s 

a 

H3H 

2.0 

6.36 

0 

A-S 

2a 


3.0 

5.11 

0 

s 

a 

b' 

3.0 

0 

1.74 

s 

2a 

2b* 



0 = 

60° 



G ;6 

0 

19.28 

s 

2a 

HI 

1.0 

0 

9.25 

s 

2a 


1.2 

16.22 

0 

A-S 

2a 


1.6 

12.32 

0 

s 

a 


2.0 

0 

4 

s 

2a 


2.0 

4.0 

3.43 

s 

2a 


2.0 

7.14 

2 

s 

2a 

K3H 

2.0 

7.97 

1 

s 

a 


2.0 

8.75 

0 

s 

a 


3.0 

6.95 

0 

A-S 

a 

B3H 

3.0 

0 

2.56 

s 

2a 



Is - symmetrical about- panel midpoiiit 
A-S - autlsymmetrlcal about panel midpoint. 
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TABLE 2.- DATA. Diib'iNIHG THE BUCKLIHG MODES IN FIGUH^ 2 AND 3 
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TABLE 3.- DATA DEJj'lWIKG THE HJCKLING MODES IR EIGUEE 5 



p/N 

q/M 

J/N 

J/M 

0 = 0°, 30° 

k 

1 

1 

0 

0 

5 

1 

1 

0 

0 

6 

1 

1 

0 

0 

0 = 45° 

k 

1 

1 

0 

0 

5 

1 

1 

0 

0 

.6 

1 

1 

0 

0 

7 

1 

1 

1 

1 

8.5 

1 

1 

3/2 

1 

0 = 60° 

4 

1 

1 

0 

0 

5 

1 

1 

1 

1 

6 

1 

1 

3/2 

3/2 

7 

1 

1 

2 

2 

8.5 

1 

1 

3 

3 




MCA TN 2392 


0 .5 

Figure 


1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

$ 


Theoretical hucklliig stressee for parallelograiii-shaped panels 
for stresB acting parallel to a set of sides. 
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UmilHH 



Results of energy solutions of higher accuracy for 
<j) of 30* 45* and 60**, respectively 



(Td) Con5)arison of variouB solutions. 

Figure 5 -- Stability of equal-sided skew panels under simultaneous action 

of stresses in x- and y-directions . 
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